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We present preliminary numerical results from a lattice study of the two-dimensional 0(3) non- linear sigma 
model. In the continuum this model possesses TV = 2 supersymmetry. The lattice formulation we use retains an 
exact (twisted) supersymmetry except for a soft breaking associated with a Wilson mass term needed to remove 
the doubles. Our numerical results show that the partition function is independent of coupling as predicted by 
supersymmetry. 



1. General 2D Sigma Model 

As was first noted by Polyakov [I], there 
exists a deeply rooted analogy between four- 
dimensional Yang-Mills theories and two- 
dimensional sigma models. Indeed the latter 
can serve as perfect theoretical laboratories to 
test methods and approaches developed for solv- 
ing the problems of these far richer and more 
complicated theories. This is particularly true 
for the supersymmetric extensions of these the- 
ories 0. It is clearly important to formulate 
such theories on the lattice and a greal deal of 
progress has been made in this area recently 
for models possessing extended supersymmetry 
|3I4I5I6| . These formulations all rely on the idea 
of twisting - a technique pioneered by Witten [7] 
in the context of topological field theory. These 
twisted formulations possess a scalar nilpotent 
supercharge Q, which with care, may often be 
preserved on the lattice. This symmetry can 
help prevent relevant susy breaking counterterms 
from appearing in the effective action and thus 
may avoid the fine tuning problems of generic 
discretizations of supersymmetric models. As 
was discussed in Simon Catterall's talk in this 
conference the twisted action is typically Q exact 
and in the case of the 2D sigma model takes the 



form 



S = aQjr hli {<!>) -\g ij B? 



where the action of the twisted charge on the 
fields is given by 



Q4> % = A 2 
QX i = 

Qmn = 



QB ip = [B m Tl k \«~- Vm R{ lk \ l \ k 



Here (j> l (a) is a set of commuting fields which act 
as coordinates on some iV-dim target space with 
metric gij{4*)- A l (cr) are scalar grassmann fields, 
Vitt( a ) are vector grassman fields and B^a) are 
another commuting vector fields. If we carry out 
the Q- variation and integrate out B we find 

S = aj^(^g l] N^N^~f ]l ^ k N^X k (1) 
+ \R j i m krf%\ m \ k ) (2) 
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To generate a kinetic term for the <j> fields it is 
necessary to choose 

iV v = <9V 

Finally to make contact with usual sigma model 
we need to impose self-duality conditions on the 



1 



2 



vector fields in the model. This amounts to re- 
placing r\ by P^rj where 



and J is a covariantly constant matrix VfeJj 
I). Manifolds possessing such a tensor field are 
termed Kahler. It can be shown that the resultant 
model can we written in complex coordinates as 

S = a J ' d 2 a {2h + -g ll d+(t) I d-(jy I 

- h+'g^D-X 7 - h + -g 1] T 1 _D + \ J 
JJ\J\J 
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2. 2D Sigma Model on the lattice 

It is trivial to see that the action is invariant 
under the twisted supersymmetry when I replace 
the continuum derivatives by symmetric finite dif- 
ferences, as the continuum Q-symmetry makes 
no reference to derivatives of the fields. How- 
ever the kernel of the (free) lattice Dirac opera- 
tor constructed this way contains both fermionic 
and bosonic doubles which have no continuum in- 
terpretation. To remove these we need to add a 
Wilson term to the lattice action. Of course it is 
not obvious that the addition of such a term is 
compatible with the topological symmetry in the 
case of a curved target space. However, in [5] it 
was shown that indeed it is possible to add po- 
tential terms to the twisted sigma models while 
maintaining the topological symmetry. The pos- 
sible terms are 

AS = c 2 V T Vi + ^X'ViVjX 7 - hi+-i 1 I + V I V 1 r 1 1 

Here, V 1 is a holomorphic Killing vector and c 2 
an arbitrary parameter. A Wilson term would 
correspond to the choice V 1 = imw4' 1 where 



m w = l - (A+A_ 



A+A: 



Many Kahler manifolds possess such a holomor- 
phic Killing vector (for example the CP N theo- 
ries). Actually Q 2 is no longer zero in such models 
but yields a diffeomorphism along the Killing vec- 
tor. For a reparametrization invariant continuum 



action this will be an exact symmetry. However 
on the lattice (where exact reparametrization in- 
variance is broken) it will yield a symmetry break- 
ing term of 0(a 3 ). As discussed in jS] such a term 
will not affect the renormalization of the model 
and we expect that the model will flow to a con- 
tinuum theory with full N = 2 SUSY in the con- 
tinuum limit. 

3. Lattice 0(3) nonlinear sigma model 

The usual 0(3) supersymmctric sigma model 
is perhaps the simplest example of a model with 
Kahler target space which may be realized as a 
twisted model. In this case the metric, connection 
and curvature are easily verified to be 
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where p = 1+uu. In this case the supersymmetric 
lattice action including Wilson terms takes the 
form 



S = 



J2 
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(Aiu) + c.c]rj 
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To proceed further it is convenient to introduce 
an auxiliary field a to remove the quartic fermion 
term. Explicitly we employ the identity 



(7) 



where N is the number of lattice sites. Thus the 
partition function of the lattice model can be cast 
in the form 



Z 



DuDaDrjDipe 



-S(u,(T 7 rj 7 X) 



(8) 
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where the action is now given by 

-^A+mA^m H — 7-{mwu){m\yu) 



5 = 



(9) 



The anticommuting fields are assembled into 
Dirac spinors 



* = 



A 



and the Dirac operator M(u, a) in the chiral basis 
is 



M(u,a) 



mw 



2p 2 



-my/u + c.c 



-D 



^mwu + c.c 



where the covariant derivative is modified to in- 
clude a coupling to the auxiliary field a 



= a? - 7 ( A +") 



(10) 



To simulate this model we have to reproduce the 
fermion determinant arising after integrating out 
the anticommuting fields. This leads to an effec- 
tive action of the form (we observe the determi- 
nant to be positive definite) 

S = aS B (u, a ) - ^Tr In (a 2 A^ (u, a)M(u, a)) (1 1) 

where Sb{u,cj) denotes the local bosonic pieces 
of the action and we have shown the depen- 
dence on coupling a explicitly. We have employed 
a Langevin algorithm with a stochastic fermion 
force estimator to represent this determinant. In 
this case physical quantities exhibit systematic er- 
rors of order the step size 0(t). 

Using the invariance of the action under the 
twisted Q, one can easily show that the partition 
function should be independent of the coupling 
constant a. A stringent test of this a indepen- 
dence of Z is gotten by measuring the expectation 
value of Sb ■ It should be clear that 
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Figure 1. a < S B > /2V 



where V denotes the number of lattice sites. The 
results are given in figH which shows data for 
^7 < Sb > from runs at a variety of couplings a 
for a lattice of size 12x12 for different time steps 
t. This is in a good agreement with the theoret- 
ical prediction above. Notice that the small de- 
viations from unity at small a's are due to small 
Q-symmetry breaking associated with the Wil- 
son term, whereas those for large a's are due to 
the systematic errors (0(t)) in the Langevin al- 
gorithm. Further numerical studies of Ward iden- 
tities and other important quantities are under- 
way. They should give us a better understanding 
of the theory and the possibility of breaking of 
the supersymmetry. 
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